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This is the users manual of the black- hole event generator BlackMax [T], which simulates the ex- 
perimental signatures of microscopic and Planckian black-hole production and evolution at proton- 
proton, proton-antiproton and electron-positron colliders in the context of brane world models with 
low-scale quantum gravity. The generator is based on phenomenologically realistic models free of 
serious problems that plague low-scale gravity. It includes all of the black-hole gray-body factors 
known to date and incorporates the effects of black-hole rotation, splitting between the fermions, 
non-zero brane tension and black-hole recoil due to Hawking radiation (although not all simultane- 
ously). The main code can be downloaded from [2J. 

1 Introduction and New Features 

BlackMax is a very versatile semi classical and quantum black-hole generator which simulates a 
number of different extra dimension models and black-hole evolution scenarios. It also gives the 
user the possibility to set many of the parameters which influence the formation and decay of black 
holes. The manual intends to explain the various parameters the user can set. 

After the release of BlackMax version 1.01 we received comments, bug-reports and requests for new 
features from the user's community, which we tried to incorparate as much as possible. The new 
features of BlackMax version 2.00 are listed below and will be discussed in the upcoming sections in 
more detail. 

1. BlackMax is able to simulate black-hole production in proton-proton, protron-anti-proton and 
electron-positron collisions. 

2. The user is able to choose among different Planck scale conventions or use his/her own con- 
vention. 

3. Graviton emission can be simulated. 

4. The user can define if baryon, lepton numbers or flavour should be conserved or not. In case 
the user does not conserve baryon numbers BlackMax events cannot be hadronised 
with Pythia or HERWIG. 

5. The parton density functions of the LHAPDF package [3] are now included into BlackMax. 

6. The user can choose to use the Yoshino-Rychkov cross section enhancement factors and sim- 
ulate the energy loss before the formation of the event horizon as described in ^l- These 
corrections to the cross section and the energy loss are only applicable for non-tension models 
and BlackMax will not apply any corrections or modifications to the energy loss if the user 
turns the Yoshino-Rychkov corrections on for a positive tension brane scenario. 
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7. The default BlackMax output is in the LHA format [S]. 

8. Graviton and photon emission is now also included into the balding phase before the black-hole 
is formed. If the user has turned on graviton emission and choses to use the Yoshino-Rychkov 
suppression factors or sets the mass loss factor to a non-zero value, BlackMax sheds this 
energy/mass by emitting two gravitons before the formation of the black-hole. If the user has 
disabled graviton emission BlackMax will shed this energy by emitting two photons. 

9. If the user simulates the split-fermion scenario BlackMax takes into account the position of 
the black-hole remnant when it calculates which particles it should decay into during the final 
burst phase. 

10. Users should note that the output format has changed with respect to the description in [T] 
and is described in section 18.21 

2 Installation 

The most up-to-date source code and TarBall can be downloaded from: 
http : //pro j ects . hepf orge . org/blackineLx/| 

Having downloaded the zipped tar file you will need to unzip it, then extract the files and make the 
executable: 

gunzip BlackMax-2 . 00 . tar . gz 
tar -xvf BlackMax-2. 00. tar 

Before compilation you will need to check the compiler version of gcc you are using. This is because 
the latest gcc compiler version (4.1.2) has changed the names of some system libraries needed to 
compile Fortran with C code. The download is configured to use gcc version 4. If you have an older 
gcc version e.g. 3.4.6 then you will need to modify the BlackMax Makefile. You can check your 
compiler version by doing the following: 

gcc — version 

Which will generate output like this: 

gcc (GCC) 3.4.6 20060404 (Red Hat 3.4.6-10) 
Copyright (C) 2006 Free Software Foundation, Inc. 

You will need to change the Fortran system library names in the Makefile in case your compiler 
version is an older one. Do this by uncommenting the following lines in the Makefile 

F77LIB =g2c 
F77C0MP=g77 

You are now ready to compile BlackMax. 
There are three different ways to run BlackMax: 

1. in standalone mode - no additional libraries required 

2. accessing PDFs from LHAPDF 

3. accessing PDFs from LHAPDF and simultaneous hadronisation from Pythia 
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To produce the exectuable for each method simply requires a diflFerent compilation/linking step and 
is described below. In all three options the default format of the event output is the Les Houches 
Accord format ^5j. This text file can be used as input to HERWIG/Pythia to hadronise the BlackMax 
events at a later datfjil 

2.1 To Run in Standalone Mode 

In this version the proton parton densities are taken from CTEQ6m which are packaged with Black- 
Max. After unpacking simply do: 

gmake BlackMaxOnly 

modify parameter.txt to pick one of the 41 CTEQ6m PDF sets that has been bundled with Black- 
Max, for example: 

choose_a_pdf _f ile(200_to_240_cteq6)0r_ >10000_f or_LHAPDF 
200 

Then run the executable 
BlackMax > &! out 



2.2 To Run with LHAPDF 

This version of the executable uses the proton parton densities from the LHAPDF library which you 
will need to download and install from here: 



http : //pro j ects . hepf orge . org/lhapdf /| 



For LHAPDF please ensure you install the package in a directory where you have write permission. 
You can do this by specifying an installation directory - for more info see the LHAPDF manual. 
Edit the BlackMax Makefile and insert the library locations. Ensure that your LD_L1BRARY_PATH 
environment variable includes the location of your newly built LHAPDF library: 

export LD_L1BRARY_PATH = $LD_L1BRARY_PATH : /data/rizvi/atlas/lhapdf - 5.3.0/lhapdf /lib 
export LHAPATH = /data/rizvi/atlas/lhapdf - 5.3.0/lhapdf /share/lhapdf/PDFsets 

Ensure you have chosen a valid PDF set in parameter.txt, for example to choose the LHAPDF 
partons from the HI PDF2000 fit of HERA data: 

choose_a_pdf _f ile(200_to_240_cteq6)0r_ >10000_f or_LHAPDF 
70050 

After unpacking the source files do: 
gmake BlackMax 

Once the executable is created you can run the program: 
BlackMax > &! out 



^Only if baryon numbers are conserved, please see section 



bullet item 29. 
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2.3 To Run with Simultaneous Pythia Hadronisation 



In order the hadronise the events during the generation job BlackMax comes with an interface to 
Pythia. The user can only run BlackMax simultaneously with Pythia if s/he is conserving baryon 



numbers (see section 8.1 bullet item 29). In order to generate fully hadronised events you will need 



to download and install the latest versions of LHAPDF and PYTHIA. They are available at: 



http://www 


hepf orge 


org/downloads/pythia6 


http://www 


hepf orge 


org/downloads/lhapdf 



BlackMax has been tested wth Pythia 6.4.10 and LHAPDF 5.3.0. Install Pythia and LHAPDF 
according to the instructions. With Pythia, create the libraries, but remove the following dummy 
routines from Pythia: 

upinit.f 
upevnt.f 
pdfset.f 
structm.f 



You will also need to remove the mention of the pdfset.f routine from the Pythia Makefile. The 
four routines above are all dummy routines which actually exist in LHAPDF. Edit the BlackMax 
Makefile and insert the library locations. Ensure that your LD_LIBRARY_PATH environment variable 
includes the location of your newly built Pythia and LHAPDF libraries: 

export LD_LIBRARY_PATH = $LD_LIBRARY_PATH : /data/rizvi/atlas/lhapdf - 5.3.0/lhapdf /lib 
export LHAPATH = /data/rizvi/atlas/lhapdf - 5.3.0/lhapdf /share/lhapdf/PDFsets 

Then create the BlackMax executable using the target "all" which will link to the Pythia and 
LHAPDF hbraries: 

gmake all 

Ensure you have chosen a valid PDF set in parameter.txt, for example: 

choose_a_pdf _f ile(200_to_240_cteq6)Dr_ > lOOOO.f or.LHAPDF 
10050 

then run the exectuable: 
BlackMax > &! out 



3 Black-Hole Production 

We assume that the fundamental quantum-gravity energy scale is not too far above the elec- 
troweak scale. Consider two particles colliding with a center-of-mass energy Eqm- They will also 
have an angular momentum J in their center-of-mass (CM) frame. By the hoop conjecture 0, if the 
impact parameter, h, between the two colliding particles is smaller than the diameter of the horizon 
of a (d -I- l)-dimensional black-hole (where d is the total number of space-like dimensions) of mass 
M — EcM and angular momentum J, 

h<2rh{d,M,J), (1) 



4 



then a black-hole with will form. The cross section for this process is approximately equal to the 
interaction area 7r(2rft)^. 

In Boyer-Lindquist coordinates, the metric for a ((i-l-l)-dimensional rotating black- hole (with angular 
momentum parallel to the uj in the rest frame of the black-hole) is: 



4~d 



sm f I r + a 
2asin2 0-'"^'"'' 



.2 n 



i-d 



dtd4> - ^^^dr' 



E(r,6l) " A 
- 9)d9^ - r^cos^ed'^-^n 

where is a parameter related to mass of the black-hole, while 

T, = r^ + a^cos^e 

and 

A = + - iir^-''. 

The mass of the black-hole is 



(2) 

(3) 
(4) 

(5) 

(6) 

r(d/2) 

is the hyper-surface area of a (d — 1) -dimensional unit sphere. The higher-dimensional gravitational 
constant is defined as 



and 



is its angular momentum. Here, 



J 



Ad-i 



2Ma 

d-l 



Gd 



(8) 



(The user is able to choose other Planck scale conventions than the one in equation [s] see Section |8| 
The horizon occurs when A = 0. That is at a radius given implicitly by 



M) - 



Here 



^id) ^ ^l/(d-2) 



(9) 



(10) 



is the Schwarzschild radius of a (d -I- l)-dimensional black-hole, i.e. the horizon radius of a non- 
rotating black-hole. Equation [T0| can be rewritten as: 



where 



k{d) 



2<i-3^(d-6)/2J 



d-l 



The Hawking temperature of a black-hole is 



T, 



H 



d~2 



(11) 

(12) 
(13) 
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If two highly relativistic particles collide with center-of-mass energy Ecu, and impact parameter 
h, then their angular momentum in the center-of-mass frame before the collision is Lin — bEcul'^- 
Suppose for now that the black-hole that is formed retains all this energy and angular momentum. 
Then the mass and angular momentum of the black- hole will be Mj„ = -Eqm and Ji„ = L^^- A 
black-hole will form if: 

h < 5„ax = 2r[^'> {EcM.braaxEcM/2) . (14) 

We see that b^ax is a function of both Ecu and the number of extra dimensions. 



We can rewrite condition ( 14 1 as 



There is one exception to this condition. In the case where we are including the effects of the brane 
tension, the metric (and hence gray-body factors) for a rotating black-hole are not known. In this 
case we consider only non-rotating black-holes. For the model with with non-zero tension brane, the 
radius of the black-hole is defined as 

^h=^, (16) 
with B the deficit-angle parameter which is inverse proportional to the tension of the brane B = 



1 — (see equation 8 in [5]). Therefore, for branes with tension 



bTar"{EcM,d) = 2ri'\EcM). (17) 
Also, for branes with positive tension only the d = 5 metric is known. 



4 Simulated Scenarios in BlackMax 

BlackMax is able to simulate several different extra dimension models, black-hole scenarios. 

• Non-rotating black-hole on a tensionless brane: For a non-rotating black-hole, we used previ- 
ously known gray-body factors for spin 0, | and 1 fields in the brane, and for spin 2 fields {i.e. 
gravitons) in the bulk. 

• Rotating black-hole on a tensionless brane: For rotating black-holes, we used known gray-body 
factors for spin 0, ^ and 1 fields on the brane. The correct emission spectrum for spin 2 bulk 
fields is not yet known for rotating black-holes, we currently do not allow for the emission of 
bulk gravitons from rotating black- holes. 

• Non-rotating black-holes on a tensionless brane with fermion brane splitting: In the split- 
fermion models, gauge fields can propagate through the bulk as well as on the brane, so we 
have calculated gray-body factors for spin and 1 fields propagating through the bulk, but 
only for a non- rotating black- hole for the split-fermion model [9]. 

• Non-rotating black-holes on a non-zero tension brane: The bulk gray-body factors for a brane 
with non-zero tension are affected by non-zero tension because of the modified bulk geometry 
(deficit angle). We have calculated gray-body factors for spin 0,1 and 2 fields propagating 
through the bulk, again only for the non-rotating black-hole for a brane with non-zero tension 
and d = 5. 

• Two particle final states: We use the same gray-body factors as a non-rotating black-hole to 
calculate the cross section of two-particle final states (excluding gravitons according to [10^). 
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5 The Black-Hole Formation 



The formation of the black- hole is a very non- linear and complicated process. BlackMax gives the 
user the possibility to set parameters which effect the formation phase of the black-hole by hand or 
use corrections calculated by Yoshino and Rychkov ^ . 

If the user would like to set the formation parameters her/himself BlackMax assumes that, before 
settling down to a stationary phase, a black-hole loses some fraction of its energy, linear and angular 
momentum. These losses are parameterize by three parameters: 1 — /e, 1 — fp and 1 — /l. Thus 
the black-hole initial state that we actually evolve is characterized by 

E ~ EinfE', 
Pz — Pzinfp', (18) 

J' ~ Linfi^; 

where Ein, Pzin and Lin are initial energy, momentum and angular momentum of colliding partons, 
while fs, fp and /l are the fractions of the initial energy, momentum and angular momentum that 
are retained by the stationary black-hole. 

6 Black-Hole Evolution in BlackMax 

The Hawking radiation spectra are calculated for the black-hole at rest in the center-of-mass frame 
of the colliding partons. The spectra are then transformed to the laboratory frame as needed. In all 
cases we have not (yet) taken the charge of the black-hole into account in calculating the emission 
spectrum, but have included phenomenological factors to account for it as explained below. 

6.1 Electric and Color Charge Suppression 

A charged and highly rotating black-hole will tend to shed its charge and angular momentum. Thus, 
emission of particles with charges of the same sign as that of the black-hole and angular momentum 
parallel to the black-hole's will be preferred. Emission of particles that increase the black-hole's 
charge or angular momentum should be suppressed. The precise calculation of these effects has not 
as yet been accomplished. Therefore, to account for these effects we allow optional phenomenological 
suppression factors for both charge and angular momentum. 

The following charge-suppression factors can currently be used by setting parameter 
charge_suppression (cf. section |8]) equal to 2. 

= exp(CQg'"'g^™) (19) 
F! = exp(C3cfcr) a = r,b,g. (20) 

Q''^ is the electromagnetic charge of the black-hole, Q*^™ is the charge of the emitted particle; c^, 
is the color value for the color a, with a = r,b,g, of the black-hole, and cj^™, is the color value for 
the color a, with a = r,b,g, of the emitted particle, and Cs are phenomenological suppression 
parameters that are set as input parameters of the generator. 

We estimate — 0{aem) and Ca — 0{as), where aem and as are the values of the electromag- 
netic and strong couplings at the Hawking temperature of the black-hole. Note that we currently 
neglect the possible restoration of the electroweak symmetry in the vicinity of the black-hole when 
its Hawking temperature is above the electroweak scale. Clearly, since aem — 10~^ we do not ex- 
pect electromagnetic (or more correctly) electroweak charge suppression to be a significant effect. 
However, since as(lTeV) ~ 0.1, color suppression may well play a role in the evolution of the 
black-hole. 
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6.2 Angular Momentum Suppression 

Since the TeV black-holes are quantum black-holes, the gray-body factors should really depend on 
both the initial and final black-hole parameters. The calculation of the gray-body spectra on a fixed 
background can cause some problems. In particular, in the current case, the angular momentum of 
the emitted particle (as indeed the energy) may well be comparable to that of the black-hole itself. 
There should be a suppression of particle emission processes in which the black-hole final state is very 
different from the initial state. We therefore introduce a new phcnomcnological suppression factor, 
parameter L_suppression, to reduce the probability of emission in which the angular momentum 
of the black- hole changes by a large amount. 

If parameter L_suppression is equal to 1 (cf. section [s]), BlackMax does suppress the increase 
the angular momentum of the black- holes. There are three angular momentum suppression models 
BlackMax is able to simulate. 

1. AArea suppression: 

If the user sets L_suppression equal to 2 the AArea suppression model is used where the 
suppression factor is defined as 

= exp(a(rf (i + Mflri\tf - 1)), (21) 

here At is the next time step in the simulation. 

2. Jbh suppression: By setting L_suppression equal to 3 the Jbh suppression model is used and 
the angular suppression factor is defined as 

i^^=exp(-CL|J'^i + A<)|). (22) 

3. A Jf,/i suppression: For the l^Jhh suppression model L_suppression has to be equal to 4 and 

= exp(-CL| J'^t + At) - J'"'(i)|). (23) 

We might expect (^^ ~ 1, however there is no detailed theory to support this; as indeed there is no 
detailed theory to choose among these three phenomenological suppression scenarios. 



7 Final Burst of Black-Holes in BlackMax 

In the absence of a self-consistent theory of quantum gravity, the last stage of the evaporation cannot 
be described accurately. Once the mass of black- hole becomes close to the fundamental scale M*, 
the classical black-hole solution can certainly not be used anymore. We adopt a scenario in which 
the final stage of evaporation is a burst of particles which conserves energy, momenturrj^ and all of 
the gauge quantum numbers. For definiteness, we assume the remaining black-hole will decay into 
the lowest number of Standard-Model particles that conserve all quantum number, momentum and 
energy. 



Black-holes do not conserve global quantum numbers, like flavor, baryon or lepton numbers [TT]. If 
the user would like to conserve global quantum numbers s/he can set this in the parameter-flle (see 
section 8.11. 



8 Input and Output 
8.1 Input 

The input parameters for the generator are read from the file parameter.txt, see Figjl] In the 
following buUeted list an explanation is given for each of the input parameters in parameter.txt. 

^Black holes are not elementary particles and do not conserve spin and angular momentum separatly. Hence spin 
is not separately conserved in BlackMax, only the total angular momentum is. 
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Number_o f_sirtiulations 
10 

lncoming_part icle { 1 : pp_2 : ppbar_3 : ee+) 
2 

Center_of_iriass_energy_o f_incoming__particle 

14000 

M_ph{GeV) 

1000. 

def inition_of_M__pl : ( 1 :M_D_2 :M_p_3 :M_DL_4 :put_in_bY_hand) 

1 

if_def inition==4 
1 . 

Choose_a_case : 

( 1 : tensionless_nonrot ating_2 : tension_nonrot ating„3 : rot ating_nonsplit_ 

4 : Lisa_two_part icles_f inal_states ) 

1 

number_o f„extra_dimensions 
4 

number_o f_splltting_dimenslons 
1 

extradirrLension_size (1/Mpl) 
10. 

tension (parameter_of_def icit_angle : l_to_0 } 
1 . 

choose_a__pdf_file (200_to_24 0_cteq6)Or_>100_for_LHAPDF 
10050 

Minimum_mass (GeV) 
5000. 

Maxrtium_mass (GeV) 
15000. 

f ix_tirtie„step (1 : f ix_2 : no) 
2 

time_step (1/GeV) 
l.e-5 

do_yo_shino {1 :do_0 :no) 


Mass_loss_f actor (O-l . 0) 
0.00 

momentum_loss_f actor (0~1 . ) 
0.0 

Angular_momentum_loss_f actor (0~1 . 0) 
0.0 

turn_on_graviton (0 :of f_l :on) 

1 

Seed 

123589341 

V:ri-e_LHA_Output_Record?_0=NO l=Yes 

L_suppression { 1 ; none„2 ; delta_area_3 : anular_morrLentum_4 : delta_angular_momentum) 
1 

angular_momentum_suppression_f actor 
1 

charge_suppression ( 1 : none_2 :do) 
1 

charge_suppression_f actor 
1 

color_suppression_factor 
20 

split_f ermion_width ( 1 /Mpl) _and_location (f rorti-15tol5) (up_to_9extradlmensions ) 

u_quark_Right (Note:do_not_lnsert_blank_spaces) 

1.0 

10.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0 
u_quark_Left (Note : do_not_insert_blank_spaces ) 
1.0 



numl:)er_o f_conservat ion 

2 



, s 


b 


u 


c 


t 


e 




1, tau, r 


iu_e, nu_ 


, 1 


1 


1 


1 


1 








0,0,0 


0, 


. 














1 


1 


1,1,1 


1 


, 




















0,0,0 





, 




















0,0,0 





, 




















0, 0, 





, 




















0, 0, 






Figure 1: parameter.txt is the input file containing the parameters that the user can set. 
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1. Number_of .simulations: sets the total number of black-hole events to be simulated; 

2. incoming.particle (1 :pp_2 :ppbar_3 : ee+) : sets the type of incoming collision particles: 

• 1 = proton-proton collisions, 

• 2 = proton-antiproton collisions, 

• 3 = electron-positron collisions. 

3. Center_ofjnass_energy_of .protons: sets the center-of-mass energy of the colliding protons 
in GeV; 

4. M_ph: sets the value of the fundamental quantum-gravity scale (M*) in GeV; 

5. def inition_of Jt_pl : (l:M_D_2:M_p_3:M_DL_4:put_in_by_hand): this sets the definition of the 
used fundamental quantum- gravity scale; 

• 1 = Convention according to and first reference of [I3j (PDG definition). 

• 2 — Convention according |14j . which is useful in quoting experimental bounds (Giddings 
and Thomas). 

• 3 = Convention given by Dimopoulos and Landsberg; M^^^ = l/Gjj, with D the total 
number of dimensions and Gu the D-dimensional Newton gravity constant. 

• 4 = User defined convention. In this case BlackMax will read the value of the next 
parameter if _def inition==4 to set the scale. The value of if _def inition==4 is kid) in 



equation 11 Note that k{d) is depending on the number of total space dimensions and the 
user has to make sure to that the value of this parameter is appropiately updated when 
the user simulates black holes in scenarios with different number of extra dimensions. 

6. Clioose_a_case: defines the extra dimension model to be simulated: 

• 1 = non-rotating black-holes on a tensionless brane with possibility of fermion splitting, 

• 2 = non-rotating black-holes on a brane with non-zero positive tension, 

• 3 = rotating black-holes on a tensionless brane with d=5, 

• A — two-particle final-state scenario; 

7. number_of _extra_dimensions: sets the number of extra dimensions; this must equal 2 for 
branes with non-zero positive tension (Choose_a_case=2); 

8. niimber_of _splitting_dimensions: 

sets the number of extra split-fermion dimensions (Choose_a_case=l); 

9. extradimension_size: sets the size of the mini-bulli[^ in units of 1/TeV (Choose_a_case=l); 
Please also refer to the discussion under item number 28. 

10. tension: sets the deficit-angle parameter B [SJ[T5]; typical values for branes with tension are 
from 1 to 0.9. (Choose_a_case=2); 

11. choose_a_pdf _f ile(200_to_240_cteq6)0r_ >100_f orXHAPDF: The user can choose to run 
BlackMax with the bundled CTEQG.lm PDFs [16l [17] which come with the BlackMax re- 
lease (see subdirectory cteq_pdf ) or to run BlackMax with the PDFs of the LHAPDF package. 
In case the user is running BlackMax in standalone mode s/he has to set this variable in 
parameter.txt between 200 and 240 to choose one of the 41 CTEQ6.1M PDF sets. It is recom- 
mended to use 200 - the others correspond to PDF uncertainties. If the user wants to use the 
PDFs of the LHAPDF package the user must link the BlackMax executable to the LHAPDF 



library (see section 2.2). The input parameter to be given here then is the chosen LHAPDF 



ID number. The definition of the LHAPDF ID numbers can be found at [5]. 



^This is the distance between fermion branes where only gauge bosons and Higgs field can propagate in split-fermion 
brane scenario. 
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12. Minimumjnass: sets the minimum mass Mmin in GeV of the initial black-holes; 

13. Maxmumjnass: sets the maximum mass M„iax in GeV of the initial black-holes; 

14. f ix_time_step: 

• In case the user has choosen to simulate black-holes in a split fermion scenario BlackMax 
will simulate the production and evaporation of the black-holes in time steps. This 
variable influence the way the size of the time steps are calculate in the split fermion 
scenario. If set equal to 1, the code uses the parameter time_step to determine the time 
interval between events; if set equal to 2 then code tries to optimize the time step, keeping 
the probability of emitting a particle in any given time step below 10%. 

• In case the user has choosen to simulate black-holes in a non-split fermion scenario 
BlackMax will simulate the production and evaporation of black-holes in fix time steps 
if its set to equal 1 using the input of time_step to set the size of the time steps. If the 
user would like to save computation time s/he can turns this off by setting this parameter 
to equal 2. In this case the code will not calculate the location of the black-hole for each 
time step, which speeds up the event generation. 

15. time_step: defines the time interval Ai in GeV^^ which the generator will use for the black- 
hole evolution if f ix_time_step = 1; 

16. do_yo_shino(l : do_0 :no) : If the user would like to include the Yoshino-Rychkov cross sec- 
tion enhancement factors and the energy loss before the event horizon fomation j4] this pa- 
rameter needs to be set equal 1 otherwise it should be set to 0. In case it is set to 1 the 
Mass_loss_f actor, momentum_loss_f actor and Angular_niomentum_loss_f actor will be ig- 
nored by BlackMax. 

17. Mass_loss_f actor: sets the energy loss factor > > 1 as defined in equationjlSj Recom- 
mended values are between 10% to 15%. This depends also on the mimimum black-hole mass 
the user wished to simulate. If the energy loss factor is too big the probability to create a 
black-hole will be to low and BlackMax would stop the generation. 

18. inomentum_loss_f actor: defines the loss factor < /p < 1 for the momentum of initial black- 
holes as defined in equation 18 Recommended values are 10% to 15%. 



19. Angular_momentum_loss_f actor: sets the loss factor < /l < 1 for the angular momentum 



of initial black-holes s defined in equation 18 Recommended values are 10% to 15% 



20. turn_on_graviton(0 : of f _1 : on) : If this parameter is set to equal 1 BlackMax produces gravi- 
tons. Otherwise, BlackMax does not produce any gravitons. 

21. Seed: This sets the seed for the random-number generator (9 digit positive integer); 

22. Write_LHA_Output_Record(0=ND , l=Yes, 2=more detailed output): If the user sets this 
variable to 1 BlackMaxLHA.txt is written containing the full event info and output.txt will 
contain input parameters and cross section information only. If the user sets this variable 
equal to 2 BlackMaxLHA.txt is written containing full event information and output.txt will 
contain information about the chosen input parameters, the cross section, the particles which 



were emitted and the black-holes produced (see section 8.2 1 



23. L_suppression: This chooses the model for suppressing the accumulation of large black-hole 
angular momenta during the evolution phase of the black-holes (cf. discussion surrounding 
equations 21p3 1; 



• 1 — no suppression; 

• 2 = A Area suppression; 

• 3 = Jbh suppression; 
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• 4 = A J suppression; 



24. aiigular_momentum_suppression_f actor: defines the phenonienological angular-momentum 
suppression factor, (cf. discussion surrounding equation 2l|23 ); Recommended value is 0.2. 



25. charge_suppression: turns the suppression of accumulation of large black-hole electromag- 
netic and color charge during the black-hole evolution process on or off (cf. dicussion sur- 



rounding equation 19) 



• 1 = charge suppression turned off; 

• 2 — charge suppression turned on; 

26. cliarge_suppression_f actor: sets the electromagnetic charge suppression factor, C,q, in equa- 



tion 19 Recommended values is 0.2. 



27. color_suppression_f actor: sets the color charge suppression factor, ^3 in equation 19 Rec- 
ommended values is 0.2. 

28. split_f ermion_width(l/Mpl)_aiid_location(from-15tol5) (up_to_9extradimensions) In the 

first line after the name of the fermion, e.g. u_quark_Right, the user can set the width of 
fermion wave functions (in M^^ units). In the next line the user can set the centers of fermion 
wave functions (in units) in split-brane models, represented as 9-dimensional vectors (for 
non-split models, set all entries to 0). 

To satisfy the suppression of proton decay, the quark brane and lepton brane must be separated 
by at least \QW (with W the width of fermion wavefunction in extra dimension) for one extra 
dimension. For higher extra dimensions, the separation can be smaller. It can be less than 
lOW/n with n the number of extra dimension. 

To suppress n-nbar oscillation, the quarks have to be separated. The recommended value needs 
to be at least 3W/n. to bW/n, with n the number of extra dimensions. 

There is no way to satisfy all the constraints at the same time for the split-fermion case. But 
there are several papers in which the possible position for each fermion from the CP violation 
constraint are calculated. [TB] H] which we would like to refer the user to. 

29. number_of .conservation: During the balding and evaporation phase, the black hole will emit 
a certain number of quarks and leptons and acquire global quantum numbers [TIJ These 
quantum numbers are recorded by BlackMax as it reaches the final burst step, at which point 
it will generate further particle emissions so that the whole process, from beginning to end, 
will obey certain user-specified global fermion number conservation rules. 

The number_of _conservation parameter specifies the number of conservation rules which are 
to follow the subsequent d,s,b,u,c,t,e,mu,tau,nu_e,nujnu,nu_tau line in the parameter 
file. The individual conservation rules are specified as an ordered set of integral coefficients 
{df} (conservation matrix) such that the quantity '^fO.fNf is a constant, where Nf is the 
number of particular fermion flavor. 

For example, baryon number B can be conserved in BlackMax by making sure that in the flnal 
burst stage, the sum of the number of emitted u and c?-type quarks compensate for the number 
u and d-type quarks emitted in previous stages. The conservation rule is thus 

iV„ + Nd + N, + Ns + Nt + Nb = m (24) 

(recalling that quarks have a baryon number of 1/3). If this is the only additional rule, then 
the following lines would be entered into the parameter flle: 



It will also acquire angular momentum and electric and color charges, but this is not relevant for this discussion. 
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number_of _conservation 
1 

d,s,b,u,c,t,e,mu,tau, nu_e , nu_mu , nu_tau 
1,1,1,1,1,1,0,0,0,0,0,0 

The zeroes for the leptons indicate that they play no role in baryon number conservation. If, on 
the other hand, the user wishes to conserve lepton number L, the quark and lepton coefficients 
are swapped: 

number_of _conservation 
1 

d,s,b,u,c,t,e,mu,tau, nu_e , nu_mu , nu_t au 
0,0,0,0,0,0,1,1,1,1,1,1 

The two rules can be combined in the parameter file to conserve B and L separately: 

number_of _conservation 
2 

d,s,b,u,c,t,e,mu,tau, nu_e , nu_mu , iiu_tau 

1,1,1,1,1,1,0,0,0,0,0,0 

0,0,0,0,0,0,1,1,1,1,1,1 

Conservation of B — L takes the following form: 

^(7V„ + Nd + N, + Ns + Nt + Nb) ~ (N, + N,^ + iV^ + 7V^,, +Nr + N,J = B-L (25) 

from which parameters follow: 

number_of _conservation 
1 

d,s,b,u,c,t,e,mu,tau, nu_e , nu_niu , iiu_tau 
1,1,1,1,1,1,-3,-3,-3,-3,-3,-3 

Finally, if all fermion flavors are to be conserved individually, twelve equations need to be 
specified: 

nmnber_of _conservation 
12 

d,s,b,u,c,t,e,mu,tau, nu_e , nu_mu , nu_tau 

1,0,0,0,0,0,0,0,0,0,0,0 

0,1,0,0,0,0,0,0,0,0,0,0 

0,0,1,0,0,0,0,0,0,0,0,0 

0,0,0,1,0,0,0,0,0,0,0,0 

0,0,0,0,1,0,0,0,0,0,0,0 

0,0,0,0,0,1,0,0,0,0,0,0 

0,0,0,0,0,0,1,0,0,0,0,0 

0,0,0,0,0,0,0,1,0,0,0,0 

0,0,0,0,0,0,0,0,1,0,0,0 

0,0,0,0,0,0,0,0,0,1,0,0 

0,0,0,0,0,0,0,0,0,0,1,0 

0,0,0,0,0,0,0,0,0,0,0,1 

It should be noted that not all the possibilities the user can enter as rules are realized in 
BlackMax. If the user specifies a rule which BlackMax cannot implement, it will print an error 
message and halt. 

Figure [l] shows an example where number_of _conservatioii is set to 2 and there are 6 lines 
in the conservation matrix. Only the first two lines of the conservation matrix in Figure [T] will 
be read and the rest will be ignored. 



13 



8.2 Output 



The BlackMax code produces 3 types of output: basic information printed to screen, the Black- 
MaxLHArecord.txt file, and the output.txt file shown in Figs. 2, 3, and 4. This is controlled by the 
parameter Write_LHA_Output_Record which can be set to 0, 1, or 2 for increasingly detailed output. 

The screen dump provides basic version and input settings, cross section and timing information 
and can be redirected to a file. This information is always given. For parameter value output.txt 
is written and includes a complete list of all input parameters and the calculated cross section only. 
For parameter value 1 (default) the BlackMaxLHArecord.txt file is additionally created containing 
all input parameters, followed by an event- by-event record in LHA format [5]. This file may be used 
as input to any other LHA compliant MC e.g. for further hadronisation by Pythia or HERWIG. 
Finally for parameter value 2 the output.txt file then also contains much more detailed event-by- 
event information for each emission step in the black-hole decay. The information is tagged by an 
ID word (Begin,Parent, Pbh, trace, Pem, Pemc or Blast). Note that the energy in the output.txt file 
is the bulk energy of the particles where as the energy in BlackmaxLHArecord.txt is the observed 
energy as defined in equation 57 in [1^. 

• Parent: identifies the partons whose collision resulted in the formation of the initial black-hole 
(see Fig. [5]). 

— column 1: identifies the black-hole; 

— column 2: PDGID code of the parton; 

— column 3: electric charge of parent parton in 3Q; 

— column 4-6: color-charge vector components of the parent parton; 

— column 7: energy of the parton in GeV; 

— columns 8-10: brane momenta of the parton in GeV; 

— colunm:ll: mass of the parton in GeV. 

• Pbh: contains the evolution of the charge, color, momentum and energy of the black-holes, 
and, for rotating black- holes, their angular momentum (cf. Fig. |6|. 

— column 1: identifies the black-hole; 

— column 2: time at which the black-hole emitted a particle; 

— column 3: PDGID code of a black-hole; 

— column 4: three times the electromagnetic charge of the black-hole; 

— columns 5 to 7: color-charge vector components of the black-hole; 

— columns 8: energy of the black-hole in the laboratory frame; 

— columns 9 to 11: brane components of the black-hole momentum in the laboratory frame; 

— columns 12 to (8-|-d): bulk components of the black-hole momentum; 

— column (9-|-d): mass of the black-hole; 

— column (10-|-d): angular momentum of the black-hole, in the case of rotating black-holes; 
empty otherwise. 

• trace: contains the evolution history of the black- holes' positions (cf. Fig. [t]): 

— column 1: identifies the black-hole; 

— column 2: the times at which the black-hole emitted a particle; 

— columns 3-5 are the brane components of the black-hole position vector when the black- 
hole emitted a particle; 
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Welcome to: BLACKMAX 

A TeV-Scale-Gravity Final-State Generator 

De-Chang Dai, Glenn Starkman, Dejan Stojkovic, 
Cigdem Issever, Eram Rizvi, Jeff Tseng 



De-Chang Dai 
Cigdem Issever 
Eram Rizvi 



de-chanq.dai^case. edu 

c . isseverl@physics . ox . ac . uk 

e . rizviSqmul . ac . uk 



See arXiv:0711.3012 [hep-ph] 

Release: 
$Narae: REL_1_07 $ 

$Id: BlackMax.cv 1.59 2008/11/03 15:56:58 rizvi Exp $ 
$Date: 2008/11/03 15:56:58 $ 



Initialising 

Reading Parameters 



***** Nevent s = 10 

***** Beam Particles = PP 

***** Centre of Mass Energy = 14000.000000 GeV 

***** M_ph = 1000.000000 GeV 

***** Model Type = No Tension, Non-rotating BlackHole 

***** Number Extra Dimensions = 4 

***** Number Split Dimensions = 1 

***** Size Extra Dimensions = 0.010000 

***** Tension Parameter = 1.000000 

***** Input PDFs = 200 

***** Minimum BH Mass = 5000.000000 

***** Maximum BH Mass - 15000.000000 

***** Time Step Fixed? = No 

***** Time Step - l.OOOOOOe-05 (1/GeV) 

***** Do YoShino suppression? - No 

***** Mass Loss Factor = C . IGOOOO 

***** Momentum Loss Factor - 0.000000 

***** Ang. Momentum Loss Factor = 0.000000 

***** Gravitons in FinalBurst? = Yes 

Table file = cteq_pdf /ctq61 . 00 .tbl 
InitLHA: Finished ipdf- 200 
Done initialising LHAPDF with set = 200 

***** Calculating Cross Section 

***** BlackMax cross section - 5.71466e-ll +/- 5.36096e-13 barn 
***** BlackMax will write output to LHA compatible text files 

***** Finished writing BlackMax LHA Header information - keep file open for event info. 
***** starting Event Generation 



Event 


Number 


1 


Event 


Number 


2 


Event 


Number 


3 


Event 


Number 


4 


Event 


Number 


5 


Event 


Number 


6 


Event 


Number 


7 


Event 


Number 


8 


Event 


Number 


9 


Event 


Number 


10 



Closed LHA event file 



BlackMax clock timing info: 



Initialisation time = 12 . 740000 seconds 

Event Loop time = 1.220000 seconds for 10 events 
Time per event = 0.122000 seconds 



Figure 2: This is the basic information which is printed to screen and which can be redirected to 
file. It contains information about version numbers, dates, input parameters and run time. 
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Printout of choosen input parameters: 

Nurabe r_o f _s imu 1 at i on s 
10 

incoming_particle (1 :pp_2 :ppbar_3 : ee+) 
1 

Center_of_mass_energy_of_incoming_particle 
14000 .000000 



Information about cross section: 



total cross section ■■ 
cross section error ■■ 



6.745298e-08 (GeV" 
6.402274e-10 (GeV^ 



.623921e-ll b 
.490485e-13 b 



Information about generated black holes and emitted particles: 



Pbh 

O.OOOOOe+00 



O.OOOOOe+00 
O.OOOOOe+00 40 



-9.48767e-05 -3.13192e-04 O.OOOOOe+00 -2.37725e-04 
6.01785e+03 O.OOOOOe+00 O.OOOOOe+00 -2.79088e+03 O.OOOOOe+00 5.33156e+03 



1.61349e+03 
4.40436e+03 
-9.48767e-05 
5.7S609e+03 



O.OOOOOe+00 
O.OOOOOe+DO 
-3.13192e-04 
-l,17708e+02 

l,17708e+02 



. OOOOOe+00 
.OOOOOe + 00 
O.OOOOOe+00 
3.60077e+01 



1 . 61349e + 03 
-4 . 40436e + 03 
-2 . 37725e-04 
-2.57130e+03 



O.OOOOOe+00 5.15956e+03 



Pern 10 . OOOOOe + 5 -1 1 1 .24224e + 02 1 . 82547e+Dl -1 . 16050G + C2 4 . 00^ 65e + 01 D . DOOOOe + 00 4 . 20000e+DD 

Pemc 10 .OOOOOe+00 5 -1 1 2 . 61338e + 02 -1.04 661e+01 -2 . 28924e + 02 1 .25557e + 2 . OOOOOe+0 4 . 20000e+00 

Elast 10 O.OOOOOe+00 11-3 9.97744e+02 -2,46146e+02 8.86896e+02 -3.85124e+02 O.OOOOOe+00 5.10000e-04 

Elast 10 O.OOOOOe+00 5-1001 5.46214e+02 4.45826e+02 8.53813e+01 -3.03778e+02 O.OOOOOe+00 4.20000e+00 



Figure 3: output.txt: There are three parts to this file. The first part is a copy of parame- 
ter. txt followed by the information on the black-hole production cross section as inferred from the 
events in this generator run and is produced for Write_LHA_Output_Record equal to and 1. The 
third part includes information about the black-hole and the emitted particles and is produced for 
Write_LHA_Output_Record equal to 2. The first column identifies what type of information each row 
is supplying: lines starting with "Begin" contain information about the emitted particles before the 
formation of the black-hole; lines beginning with "parent" have information about the two incoming 
partons; lines beginning "Pbh" contain information on the energy and momenta of the produced 
black-holes; lines starting with "trace" describe the location of the black-holes; rows beginning with 
"Pem" contain information about the emitted particles in the lab frame; lines headed by "Pemc" 
have the information about the emitted particles in the center-of-mass frame; rows starting with 
"Elast" describe the final burst. 
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*** Welcome to: BLACKMAX 

*** Release tag : $Name: REL_1_07 $ 

*** Release id : $ld: BlackMax.c.v 1.59 2008/11/03 15:56:58 ■ 

*** Release date: $Date: 2008/11/03 15:56:58 $ 

*** Writing input parameter set... 

Numb e r_o f _ s imu 1 at i o n s 



incoming_particle (1 :pp_2 :ppbar_3 : ee+) 



Finished writing input parameter set 



.000DD00e+D3 
5 . 3605564e-01 
.OOOOOOOe+00 
506 

507 

12 



503 
502 


505 
506 
507 



7 .0000000e+D3 4 4 ^1 41 3 

. OODOOOOe+00 1 
1 .3994101e + D4 -1 . 000 DOOOe + 00 -1 . ODOOOODe + 00 

O.DOOOODOe+OD . DOOOOD De+00 2 . 4 8761 39e+03 
O.OOOOODDe+00 
2.1139359e+02 
2.1139359e+02 
1.8526639e+00 
6.3289592e+02 
5.0680883e+01 
1.4407840e+01 
3.313501De+02 
4.525887De+02 
2.0443502e+01 
1.0569136e+02 
3.1337342e+02 
1.0426138e+02 
4.5754642e+01 
5. 781541 8e+02 
1.8985992e+02 
6.8555863e+01 



503 -: 



O.DOOOODOe+00 
4.5124750e+01 
4.5124750e+01 
1.5431276e+02 
2.1118507e+02 
5.4745365e+01 
1 . 9032106e+01 
2311343e+02 
3.6006468e+02 
■2.0739245e+02 
4.8735597e+02 
■2.8014190e+02 
2.1565688e+02 
2.0297186e+00 
3.5278679e+02 
7.0120765e+01 
2.4851565e+01 



-3.5576648e+03 
-6.6818111e+01 
6. 6818111e+01 
1.3976068e+01 
3.8860434e+01 
1.1459160e+02 
1.3500616e+02 
-4 . 4979714e + 02 
-3.3742488e+02 
-4.3278013e+02 
-2.6101426e+01 
-6.7788477e+01 
-2.6021028e+01 
-8.5747871e+01 
3.8622763e+02 
-2.6320670e+02 
-6.9845149e+01 



2.4876139e+03 
3.5576648e+03 
2.2624799e+02 
2.2624799e+02 
1.5496010e+02 
6.6833114e+02 
1.5863487e+02 
1.3715454e+02 
6. 015729De + 02 
6. 6958096e+02 
4.8035012e+02 
4.9959765e+02 
4.2578787e+02 
2.4098353e+02 
1.2616848e+02 
7.9975757e+02 
3.3202848e+02 
1.0097453e+02 



3.1999412e-03 
3.2DD230ae-03 
3.0462581e-06 
3.0462581e-06 
1.2000000e+00 
7.0000000e-02 
8.0425000e+01 
4.2000000e+00 
3.2DD0132e-03 
7.0000000e-02 
4.2000000e+00 
4.2000000e+00 
4.2000000e+00 
4.2000000e+00 
8.0425000e+01 
1.7810000e+02 
1.2000000e+00 
7.0000000e-02 



. OOOOOOOe+00 
. OOOOOOOe+OO 
. OOOOOOOe+OO 
0. OOOOOOOe+00 
0. OOOOOOOe+00 
0. OOOOOOOe+OO 
0. OOOOOOOe+OO 
0. OOOOOOOe+OO 
. OOOOOOOe+OO 
. OOOOOOOe+OO 
. OOOOOOOe+00 
0. OOOOOOOe+OO 
0. OOOOOOOe+OO 
0. OOOOOOOe+OO 
0. OOOOOOOe+00 
0. OOOOOOOe+00 
0. OOOOOOOe+00 
0. OOOOOOOe+00 



0. OOOOOOOe + OO 
0. OOOOOOOe + OO 
0. OOOOOOOe + OO 
.OOOOOOOe+00 
0. OOOOOOOe+00 
0. OOOOOOOe+OO 
0. OOOOOOOe+00 
0. OOOOOOOe+OO 
O.OODDOOOe+OO 
O.OODOOOOe+OO 
O.OODDOOOe+00 
O.OODDOOOe+00 
O.OODDOOOe+00 
.OOOOOOOe+OO 
0. OOOOOOOe+00 
0. OOOOOOOe+00 
0. OOOOOOOe+00 
0. OOOOOOOe+00 



Figure 4: BlackMEixLHArecord.txt: This file has the output of BlackMax in LHA format. 





BH ID 


PDGID 


3Q 


Colours 






Energy 




Px 




Py 




Pz 




Mass 


Parent 


1 


4 


2 








1 


1 


61349e+03 





OOOOOe+00 





OOOOOe+00 


1 


61349e+03 


1 


20000e+00 


Parent 


1 


2 


2 





1 





4 


40436e+03 





OOOOOe+00 





OOOOOe+00 


-4 


40436e+03 


3 


20000e-03 


Parent 


2 


21 





1 


- 


1 


2 


58804e+03 





OOOOOe+00 





OOOOOe+00 


2 


S8804e+03 





OOOOOe+00 


Parent 


2 


21 





-1 


1 





2 


73896e+03 





OOOOOe+00 





OOOOOe+00 


-2 


73896e+03 





OOOOOe+00 


Parent 


3 


2 


2 








1 


5 


17716e+03 





OOOOOe+00 





OOOOOe+00 


5 


17715e+03 


3 


20000e-03 



Figure 5: Lines in the output file headed by the ID = Parent contain information about the initial 
partons which formed the black-hole. 





BH ID 




Time PDGID 


3Q 


Colors 






Energy 




Px 




Py 




Pz 




Pextra . . 




Mass 




Jbh 


Pbh 


1 





OOOOOe+00 


40 


4 





1 


1 


6 


01785e+03 





OOOOOe+OO 





OOOOOe+OO 


-2 


79088e+03 





OOOOOe+OO 


5 


33156e+03 





OOOOOe+00 


Pbh 


1 





OOOOOe+00 


40 


3 


1 


1 


1 


5 


76609e+03 


-1 


17708e+02 


3 


60077e+01 


-2 


57130e+03 





OOOOOe+OO 


5 


15956e+03 





OOOOOe+00 


Pbh 


1 





OOOOOe+00 


40 


1 





1 


1 


5 


23160e+03 


-5 


43908e+01 


1 


80610e+02 


-2 


09272e+03 





OOOOOe+OO 


4 


79110e+03 





OOOOOe+00 


Pbh 


1 





OOOOOe+OO 


40 


1 





1 


1 


4 


66423e+03 


-3 


43247e+02 


-3 


07262e+02 


-2 


07124e+03 





OOOOOe+OO 


4 


15364e+03 





OOOOOe+OO 


Pbh 


1 





OOOOOe+00 


40 


3 


1 


1 


1 


4 


40733e+03 


-2 


6S425e+02 


-1 


68412e+02 


-1 


97431e+03 





OOOOOe+OO 


3 


92763e+03 





OOOOOe+00 



Figure 6: Lines in the output file headed by the ID = Pbh contain the energies and momenta of the 
black-holes for each emission step. In case of rotating black-holes, the last column in the line is the 
angular momentum. 
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— columns 6 to (2+d): the bulk components of the black-hole position vector, when the 
black-hole emitted a particle. 

• Pem: contains a list of particles the black-holes has emitted during its evaporation phase (cf. 
Fig.js]): 

— column 1: identifies the black- hole; 

— column 2: the times at which the black-hole emitted a particle; 

— column 3: PDGID code of the emitted particle; 

— column 4: three times the charge of the emitted particle; 

— columns 5 to 7: color-vector components of the emitted particle; 

— columns 8: energy of the emitted particle in the laboratory frame in GeV; 

— columns 9 to 11: brane components of the momentum of the emitted particle, in the 
laboratory frame in GeV; 

— columns 12 to (8-|-d): bulk components of the momentum of the emitted particle in GeV; 

— column 9+d: mass of the emitted particle in GeV. 

• Pemc: contains the same information as Pem, but in the center-of-mass frame of the collision. 

• Begin: contains the same information as Pem, but for the particles which are emitted before 
the black-hole is formed. This happens if the user chooses to use the Yoshion and Rykhov 
factors or sets the mass-loss-factor to a non-zero value. 

• Elast: contains the same information as Pem for the particles emitted in the final decay burst 
of the black-hole but column 12 is the mass of the particle. The is no information on the bulk 
momentum since these particles have no bulk momentum. 

A Comparison between BlackMax and Charybdis 

In this section we compare the cross sections of BlackMax and Charybdis ^20j which is a very com- 
monly used black hole generatoi[^ There are several differences between BlackMax and Charybdis 
which needs to be considered. 

For a direct comparison one has to make sure that both generators use the same convention for 
the definition of the Planck Mass. In Charybdis the user can choose between the different con- 
ventions by setting the parameter MSSDEF. In BlackMax the user can switch between the different 
conventions with the help of the parameter def inition_of _M_pl. 

• PDG definiton, M^,: 

def inition.of _M_pl 1 corresponds to MSSDEF 3; 

^For the comparison in this section MSRT98LO was used as the input PDF. 



Time 



trace 
trace 
trace 
trace 
trace 



. ODODOe+00 
O.OOOOOe+00 
O.OOOOOe+00 
O.OOOOOe+00 
. OOOOOe+00 



. 48767e-05 
. 48767e-05 
. 48767e-05 
. 48767e-05 
. 48767e-05 



13192e-04 
13192e-04 
13192e-04 
13192e-04 
13192e-04 



O.OOOOOe+00 
O.OOOOOe+00 
O.OOOOOe+00 
O.OOOOOe+OO 
O.OOOOOe+OO 



37725e-04 
37725e-04 
37725e-04 
37725e-04 
37725e-04 



Figure 7: Lines in the output file headed by the ID = trace contain the location of the black-hole 
fo reach emission step. 
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BH ID 


Time 


PDC3ID 


3Q 




Colors 




Energy 




Px 




Py 




Pz 


Pextra... 




Mass 


Pem 


1 


. OOOOOe+00 


-5 


1 


-1 








2 


. 5n57e + 02 


1 


. 17708e+02 


-3. 


. 60077e+01 


-2 


.19573e+02 


. OOOOOe+00 


4 . 


.20000e+00 


Pem 


1 


.OOOOOe+00 


G 


2 


1 








5 


. 34489e+02 


-6 


. 33172e+01 


-1 . 


. 44602e+02 


-4 


. 78582e+02 


. OOOOOe+00 


1 . 


. 78100e+02 


Pem 


1 


.OOOOOe+00 


-16 














5 


. 67379e+02 


2 


.88857e+02 


4 . 


. 87872e+02 


-2 


.14790e+01 


. OOOOOe+00 


0. 


.OOOOOe+00 


Pem 


1 


.OOOOOe+00 


-6 - 


-2 


-1 








2 


.56892e+02 


-7 


.48229e+01 


-1 . 


.38850e+02 


-9 


. 69335e+01 


. OOOOOe+00 


1 . 


. 78100e+02 


Pem 


1 


.OOOOOe+00 


-1 


1 


-1 








8 


.26391e+02 


1 


, 19497e+02 


-1 . 


. 63989e + 02 


-8 


. 01093e+02 


. OOOOOe+00 


3. 


.20000e-03 



Figure 8: Lines in the output file headed by the ID = Pem contain the types of the emitted particles, 
their energies and momenta in the lab frame and the times of their emission. 





BH ID 




Time PDGID 


3Q 


Colors 






Energy 






Px 




Py 




Pz 




Pextra... 




Mass 


Blast 


1 





.OOOOOe+00 


2 


2 


1 








1 


. 69854e + ( 


)2 


1 


. 49973e + 02 


5 . 


.22441e+01 


-G 


. 0240Ge+01 





. OOOOOe+00 


3 


.20000e-03 


Blast 


1 





, OOOOOe+00 


4 


2 


1 








1 


. 73812e+( 


)2 


1 


. 32172e+02 


4 . 


. 97772e+01 


-1 


. 01301e+02 





.OOOOOe+00 


1 


.20000e+00 


Blast 


1 





, OOOOOe+00 


1 


-1 





1 





8 


. 97411e + ( 


)2 


1 


. 86142e+02 


4 . 


. 63555e+01 


-8 


.7GGG9e+02 





. OOOOOe+00 


3 


.20000e-03 


Blast 


1 





, OOOOOe+00 


5 


-1 








1 


2 


. 06777e + ( 


)2 


-1 


. 58071e+01 


1 . 


.36863e+02 


-1 


. 5413Ge+02 





. OOOOOe+00 


4 


.20000e+00 


Blast 


2 





.OOOOOe+00 


21 








1 ■ 


-1 


2 


. 09713e + ( 


)2 


-2 


. 09011e+02 


-1 . 


. 42072e+01 


9 


. 588G3e+00 





. OOOOOe+00 





.OOOOOe+00 



Figure 9: Lines in the output file headed by the ID = Elast contain the types, energies and momenta 
of particles of the final burst. 



• Giddings and Thomas, Mpi 

def inition_of _M_pl = 2 corresponds to MSSDEF = 1; 

• Dimopoulos and Landsberg, Mdl' 

def inition_of _M_pl = 3 corresponds to MSSDEF = 2. 

One important difference between the two generators is the definition of cross sections. Charybdis 
uses the cross section definition for a non-rotating black-hole 

<7ch = Tirf . (26) 



Here, is the Schwarzchild radius of the black-hole. 
But BlackMax uses 



Cbm = b^TTT^ (27) 

hd = ^ T- (28) 

(l + (^)2)^ 



which is the definiton of the cross section for rotating black-holes (see also equation ([T5|), where d 
is number of space dimensions. The variable bd is in general larger than one. Therefore one expects 
that the cross section of BlackMax is larger than the cross section of Charybdis. Ubm/cc/i is expected 
to be equal to 6^. In the tables [l] to [s] the value of (Jb7yi/<^ch is in general smaller than b% and in some 
cases it is even smaller than 1 . This difference comes from the fact that BlackMax assumes that the 
width of the extradimension is Af^^^, and Charybdis assumes that the width of extradimension is 
0. Because of the finite width of extradimensions, the cross section of BlackMax is further reduced 
compared to Charybdis. 



To make sure that the two generators indeed output the same cross section in the same situation, 
we manually put the size of extradimension to zero in BlackMax for the purpose of this comparison. 
The results are shown in the sixth column of the tables [T] to |3] and they agree with the expected 
value well and there is a small difference between 6^ and Ubm/'^ch of 3%. 
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Table 1: Md=1000 GeV, Mbh > 5000 GeV, and D is the total number of dimensions (space + time). 



D 


<Jch [pb] 


Cfem [pb] 


abm with Lextra = [pb] 




O'hm/o'c/i with Lextra — 




6 


75.20 ± 0.6968 


90.69 ± 0.8407 


99.70± 0.9128 


1.21 


1.32 


1.36 


7 


122.0± 1.126 


161.9± 1.502 


177.0±1.638 


1.32 


1.45 


1.48 


8 


172.6 ±1.590 


247.6 ± 2.304 


266. 2± 2.449 


1.43 


1.54 


1.59 


9 


225.7± 2.076 


352.7 ± 3.149 


369.0±3.285 


1.56 


1.63 


1.69 


10 


280.7± 2.579 


455.2± 4.182 


484.8± 4.419 


1.62 


1.72 


1.78 



Table 2: Aip=1000 GeV, Mt,h > 5000 GeV, and D is total number of dimensions (space + time). 



D 


o-ch [pb] 


O-fem [pb] 


abm with Lextra = [pb] 




fbrn/fcft, with Lextra = 




6 


119.4±1.106 


149.1 ± 1.375 


158. 3± 1.449 


1.25 


1.33 


1.36 


7 


172.6±1.593 


234.5±2.158 


250.4± 2.316 


1.36 


1.45 


1.48 


8 


227.8±2.098 


334.2±3.090 


351.2±3.231 


1.47 


1.50 


1.59 


9 


284.4±2.626 


450. 2± 4.010 


464.8±4.139 


1.58 


1.63 


1.69 


10 


342.1±3.144 


561.5 ±5.186 


591. 0± 5.387 


1.64 


1.73 


1.78 
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